c“:> 


UNCLASSIFIED 


FENSE  DOCUMENIATiON  CENTER 


ItNTIFIC  AND  TtCHNICAl  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  goveranent  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  i>er8on  or  corporation,  or  conveying  any  rlf^ts 
or  x>ermission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


CATALOGED  BY  DDC 

AS  AD  No. 


GJ_3L_4 


AERONUTRONIC  DIVISION 


Publication  No.  U-2148 


AERONUTRONIC  DIVISION 


RESEARCH  LABORATORY 


TECHNICAL  REPORT 


CHEMICAL  CORROSION  OF  ROCKET  LINER 
MATERIALS  AND  PROPELLANT 
PERFORMANCE  STUDIES 

Fourth  Quarterly  Technical  Summary  Report 


Prepared  for:  Advanced  Research  Projects  Agency 

The  Pentagon 
Washington  25,  D.  C. 

Under  Contract:  NOw  61-0905-c,  Task  E 

With  the  Bureau  of  Naval  Weapons 
ARPA  Order  No.  22-62 


Prepared  by: 


R.  C.  Oliver 
R.  W.  Baier 
D.  L.  Peters 
R.  W.  Sprague 


15  June  1963 


DDC 


U  JUN2il963 

iHtc'raEij  w  Tsiiyj 

TISIA  a 


AERONUTRONIC  DIVISION 


ABSTRACT 


Studies  on  the  theoretical  equilibria  resulting  from  Interaction 
of  candidate  refractory  nozzle  wall  materials  with  various  chemical  species 
resulting  from  combustion  of  rocket  fuels  have  continued.  Tantalum  carbide 
has  been  considered  In  detail,  with  composition  diagrams  and  saturation 
parameters  reported  for  Its  Interaction  at  1000  psla  with  AIF^,  BF2,  BF^, 

BOF,  BeF2i  CO,  CO2,  HCl,  HF,  H2,  H2O,  LIF  and  N2.  Similar  computations  are 
reported  for  tungsten  Interacting  with  BOF,  CO,  and  CO2,  these  results 
representing  a  revision  of  data  reported  In  the  previous  quarterly  report 
under  this  contract,  as  necessitated  by  a  recent  revision  In  the  thermo¬ 
dynamic  data  for  the  tungsten  oxides.  The  oxides  are  considered  briefly  as 
a  class  for  application  In  oxygen-containing  atmospheres,  with  results  given 
for  MgO,  Hf02  and  Th02  In  CO2.  The  stability  of  carbides  In  nitrogen  Is  also 
considered. 

The  possibility  of  extending  the  Isothermal  saturation  results  for 
pure  components  at  1000  psla  to  mixtures  at  other  pressures  has  been  considered 
In  some  detail,  with  results  that  look  promising  except  where  multi -atomic 
species  are  Important  which  are  not  considered  In  the  Individual  cases  (e.g., 
HCN  Is  Important  In  a  equlllbri’jm  with  graphite,  but  Is  not 

considered  with  the  Individual  pure  reactants).  A  "catalog”  of  ablation 
reactions  for  studies  reported  previously  Is  Included. 

Plans  for  future  work  are  described  briefly.  Certain  errata  for 
previous  reports  are  appended. 
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SECTION  1 

INTRODUCTION 


This  Is  the  fourth  quarterly  report  In  a  program  which  has  two 
parts,  the  first  being  the  theoretical  appraisal  of  various  refractory 
rocket  liner  materials  from  a  chemical  standpoint  and  the  second  being 
the  computation  of  theoretical  performance  of  certain  rocket  propellant 
combinations.  Work  In  this  past  quarter  has  been  concerned  only  with  the 
evaluation  of  refractory  nozzle  materials. 

Refractory  materials  under  study  Include  graphite,  tungsten,  and 

various  high  melting  carbides,  nitrides,  borides,  and  oxides.  Ihe  procedure 

In  use  for  comparing  these  materials  Is  to  compute  the  theoretical  amounts 

of  these  refractory  materials  required  to  "saturate"  a  unit  mass  of  reacting 

gas  or  solid  which  represents  a  component  of  a  rocket  combustion  stream;  the 

details  of  the  procedure  were  given  In  the  Volume  I  of  the  Second  Quarterly 

Report^.  As  explained  therein,  stability  Is  expressed  In  terms  of  both  an 

Isothermal  saturation  parameter  (the  value  "A")  and  an  adiabatic  saturation 

parameter  (the  value  "B").  The  adiabatic  saturation  parameter  Is  proportional 

under  certain  conditions  to  the  rate  of  attack  by  a  flowing  stream.  Both 

parameters  Include  effects  of  simple  vaporization,  which  factor  becomes 
2 

Important  even  In  the  absence  of  chemical  reactions  at  tenqieratures  greater 
than  about  3500°K,  especially  with  low  molecular  weight  gas  streams  and  with 
more  volatile  refractories.  The  computational  technique  Is  arbitrarily 
extended  to  temperatures  above  the  melting  point  without  modification,  a 
technique  which  Is  valid  In  terms  of  chemical  stability  but  obviously  not 
significant  In  terms  of  physical  stability. 
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Thermochemical  data  used  In  these  calculations  have  been  obtained 

3 

for  the  most  part  from  the  JANAF  compilations.  It  has  In  addition  been 

necessary  to  accumulate  thermodynamic  data  for  certain  heavy  metal  species 

not  available  In  these  tables,  and  this  has  been,  prior  to  this  quarter, 

represented  a  significant  portion  of  the  effort  In  this  program. 

Chemical  compatibility  studies  are  reported  herein  for  tantalum 

carbide  with  various  pure  gases  and  solids,  and  for  tungsten  In  oxygen- 

containing  gases  (with  the  exception  of  H2O).  In  addition  oxides  as  a  class 

are  discussed  briefly  for  application  In  oxidizing  atmospheres,  as  well  as 

nltrlde-carblde  equilibrium.  The  tungsten -oxygen  calculations  represent  a 

2 

revision  of  results  reported  previously  ,  necessitated  by  recent  revisions 
In  the  JANAF  tabulations  of  data  for  tungsten  oxides.  Equilibrium  systems 
In  which  H2W0^(g)  may  appear  have  not  been  reported,  as  revision  of  H2W0^(g) 
thermodynamic  data  Is  expected  In  the  near  future  . 

The  use  of  an  averaging  technique  to  predict  the  "A"  saturation 
parameter  for  real  combustion  gas  mixtures  Is  discussed. 


AERONUTI«ONIC  DIVISION 


SECTION  2 

CHEMICAL  CORROSION  STUDIES 


2.1  CHEMICAL  CORROSION  OF  TANTALUM  CARBIDE 

2.1.1  GENERAL 

The  chemical  stability  of  tantaliim  carbide  In  fourteen  dlff- 

1  2 

erent  atmospheres  was  calculated  as  described  In  previous  reports  ’ 
under  equilibrium  conditions  at  1000  psla  and  various  tenperstures ; 
the  atmospheres  were  formed  by  equilibration  of  the  nominal  species 
AlFj,  BF^,  BF2,  BOF,  BeF2,  CO2,  CO,  HCl,  HF,  H2,  H2O,  LIF  and  N2.  (The 
nominal  species  specifies  the  overall  stoichiometry;  the  species  Itself 
may,  however,  be  largely  dissociated  at  the  temperature  and  pressure 
under  consideration.)  The  results  are  presented  as  "A"  and  "B"  plots  in 
Figure  1;  equilibrium  compositions  are  given  In  Figures  2  through  6. 
Principal  species  contributing  to  corrosion  are  shown  In  Table  I. 

It  Is  apparent  from  Figure  1  that  in  many  cases  TaC  is  severely 
attacked  (B  values  greater  than  unity)  by  species  containing  "loosely" 
bound  oxygen,  such  as  CO2  and  H2O;  species  such  as  BF^  and  BeF2  which 
contain  "excess"  fluorine  are  less  damaging.  Exothermic  effects  are 
noted  for  the  species  HCl,  HF  and  N2,  and  endothermic  effects  with  AIF^, 

BF^,  BF2  and  BOF. 

The  species  HBO2  was  considered  but  again  gave  rather  meaningless 
results  as  no  gas  phase  was  found  to  be  present  below  about  4500-5000°K,  (as 
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FIGURE  1.  CHEMICAL  CORROSION  OF  TaC  BY  COMMON  COMBUSTION  PRODUCTS 
AT  1000  PSIA 
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FIGURE  2.  SPECIES  PRESENT  UNDER  EQUILIBRIUM  CONDITIONS  AT  10^0  PSIA 
FOR  TaC  AI«)  AIF^,  BF„  AND  rT 


FIGURE  3.  SPECIES  PRESENT  UNDER  EQUILIBRIUM  CONDITIONS  AT  1000  PSIA 
FOR  TaC  AND  BOF,  BeF-,  AND  CO 
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FIGURE  t  .  SPECIES  PRESENT  UNDER  EQUILIBRIUM  CONDITIONS  AT  1000  PSIA 
FOR  TaC  AND  CO  ,  HCl,  AND  HF 
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TABLE  I 

PRINCIPAL  SPECIES  CONTRIBUTING  TO  THE  ABLATION  OF 
TANTALUM  CARBIDE  AT  1000  PSIA 

Approximate 


Combustion  Principal  Species  at  the  Temperature 
Product  _ NoEalei  Wall _  Range  (°K)  Remarks 


!. 

TaFj,  C(c).  C2F2 

3000-5000 

No  gas  phase  below 
about  3000°K. 

•  BF2 

TaB2(c),  C(c) 

1500 

- 

TaFj,  TaF^,  TaF^, 

C2F2. 

C 

3500-5000 

BF^ 

TaFg,  C(c),  C2F2. 

TaF^, 

TaF^ 

2500-5000 

:  BOF 

Ta(c),  CO 

3500-5000 

Ta20g(c),  TaB2(c) 

2500-3000 

BeF2 

C(c),  TaFg,  TaF^, 

TaF3, 

S"2 

2000-5000 

CO 

TaO 

4500-5000 

CO2 

Ta205(c) 

2000-3500 

Ta(c),  TaO 

4000-5000 

HCl 

C(c) 

1500-3000 

TaCl^,  TaClj 

2000-5000 

Decrease  with  temp. 

HF 

C(c),  TaFj 

1500-5000 

Decrease  with  temp . 

TaF3,  TaF^,  C2H 

3000-4500 

Increase  with  temp . 

«2 

Ta(c) 

3000-5000 

H2O 

C(c).  Ta205(c),  Ta(c) 

CO 

1500»3000-5000 

1500-5000 

Solids  exist  In  three 
separate  temp,  ranges 

LIF 

TaFj,  C 

4500-5000 

^2 

CN,  TaN(c) 

3500-5000 
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2 

was  the  case  with  tungsten  )  and  this  result  was  considered  unreliable 
In  view  of  the  lengthy  extrapolation  used  in  preparing  the  thermodynamic 
tables  for  this  species  in  this  region.  It  accordingly  has  not  been  re¬ 
ported  in  detail. 

2.1.2  DISCUSSION  OF  INDIVIDUAL  CASES 
AlFj 

The  saturation  numbers  increase  over  the  temperature  range  in¬ 
vestigated,  with  a  fairly  constant  cooling  of  about  200°.  The  cooling  can 
be  attributed  to  endothermic  reactions  involving  formation  of  tantalum  and 
carbon  fluorides. 

BF2  and  BF^ 

The  general  behavior  of  these  systems  is  similar  to  the  AIF^ 
system.  Cooling  is  slightly  greater  in  the  case  of  BF^. 

BOF 

This  system  forms  solid  phases  which  have  a  vapor  pressure  less 
than  1000  psia  below  about  2500°K.  At  about  3000-3500°,  the  reduction  of 
the  carbide  to  the  metal  with  formation  of  CO  produces  a  substantial  cool¬ 
ing  effect  with  the  equilibrium  wall  temperature  over  1000°  cooler  than  the 
free  stream  temperatures. 


This  system  shows  essentially  a  constant  saturation  number  as  was 

2 

found  also  for  tungsten  with  this  gas  .  Basically  the  BeF2  is  dissociating 
to  BeF  and  F,  the  free  fluorine  then  reacting  with  the  condensed  phase. 

The  overall  reaction  is  endothermic  at  low  temperature  (640°  at  2500°)  and 
exothermic  at  high  temperature  (1100°  at  4500°).  It  should  be  mentioned 
that  these  results  are  based  on  the  JANAF  value  of  -88.34  kcal/mole  for  the 
standard  heat  of  formation  of  BeF  at29£PK;a  less  negative  value  which  now 
appears  probable  would  presumably  materially  alter  the  results. 
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CO 

Very  little  reaction  takes  place  over  the  whole  temperature 
range,  as  would  be  expected.  Some  cooling  occurs  at  very  high  temper¬ 
atures  (53°K  at  5000°K;  11°K  at  4500°K)  due  to  volatilization. 


The  corrosion  Is  essentially  due  to  oxygen  from  dissociation  of 
CO^.  At  teiq>eratures  below  3500°K,  Ta^O^/c  Is  formed;  at  temperatures 
somewhat  above  3500  k,  Ta/c  Is  formed.  The  latter  reaction  results  In  a 
larger  saturation  value  than  the  former.  The  reaction  becomes  strongly 
exothermic  at  high  temperature  (-211°  at  2500°,  +255°  at  3000°  and  1497° 
at  4500°). 

HCl 

This  gas  reacts  primarily  (and  exothermically)  with  TaC/c  to 
form  solid  carbon,  and  gaseous  tantalum  chlorides.  Saturation  numbers 
decrease  with  temperature  due  to  decreasing  stability  of  the  tantalum 
chlorides  and  Increase  In  the  fraction  of  chlorine  that  remains  In  the 
elemental  form.  The  equilibrium  wall  temperature  Is  some  100-400°  greater 
than  the  free  stream. 


HF 

This  gas  behaves  chemically  In  a  fashion  similar  to  HCl.  The 
heating  effect  Is  greater  at  low  temperatures  (nearly  1000°K)  and  less  at 
high  temperatures  (about  100°K  at  5000°) . 


H 


2 


Very  little  reaction  occurs  below  about  3000°K,  at  which  temper¬ 
ature  the  saturation  reaction  Is  slightly  exothermic  (23°K) .  The  saturation 
numbers  Increase  with  temperature  as  more  gaseous  hydrocarbons  and  elemental 
metal  are  formed;  the  net  heating  effect  decreases  gradually  until  at 
5000°K  the  net  reaction  Is  endothermic  by  1.3°K. 
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H2O 


Up  to  3000°K,  primary  reaction  product,  and  Is 

formed  nearly  quantitatively.  At  higher  temperatures  elemental  tantalum 
forms.  The  saturation  numbers  are  high  over  the  whole  range  Indicating 
that  H2O  Is  extremely  corrosive  to  TaC.  The  reactions  are  exothermic  In 
the  lower  temperature  region  and  strongly  endothermic  where 


stable  and  Ta/c  forms. 


LIF 

No  gas  phase  Is  computed  below  about  3500*^1%..  Gaseous  products 
Include  tantalum  and  carbon  fluorides  along  with  gaseous  lithium.  The 
overall  reactions  are  very  slightly  endothermic  (15°K-50°K)  over  the 
3500-5000°K  range. 


TaC  reacts  quantitatively  with  to  form  TaN  and  carbon  at  low 
temperatures  (see  Section  5).  At  3000°  the  reaction  is  exothermic,  so 
that  the  adiabatic  wall  is  some  400°K  hotter  than  the  free  stream.  The 
reaction  Is  exothermic  to  about  this  same  degree  to  4500°K;  at  5000°K 
the  reaction  Is  only  38°K  exothermic  due  to  increase  in  endothermic  vola¬ 
tilization  effects.  CN  Is  the  primary  gaseous  reaction  product. 


2.2  RECALCULATION  OF  THE  CORROSION  OF  TUNGSTEN  BY  BOF,  CO  AND  CO2 

2 

Since  publication  'of  the  Third  Quarterly  Report,  revised  or  new 

3 

data  sheets  have  been  Issued  for  the  JANAF  tables  on  the  WO^/c,  WO^Cg) , 
W^O^(g)  and  W^O^^Cs)  species.  As  substantial  changes  were  Involved,  com¬ 
putations  for  cases  where  these  oxides  might  be  significant  were  repeated. 

(The  tungsten-water  case  has  again  been  held  up  pending  anticipated  receipt 

4 

of  new  Information  on  H2W0^(g)  .  The  HBO2  case  was  not  recomputed  In  view 
of  the  rather  meaningless  results  obtained  before^,  with  no  gas  phase 
present  at  1000  psla  below  about  4500°K.) 

Results  are  given  In  Figures  7  and  8. 
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FIGURE  7.  CHEMICAL  CORROSION  OF  TUNGSTEN  BY  BOF,  CO  AND  CO2  AT 
1000  PSIA  (REVISED  CALCULATIONS) 
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FIGURE  {;.  SPECIES  PRESENT  UNDER  EQUILIBRIUM  CONDITIONS  AT  1000  PSIA 
FOR  TUNGSTEN  AND  BOF,  CO  AND  CO  (REVISED  CALCULATIONS) 
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Comparison  of  the  results  with  the  new  and  old  data  can  be  made 
by  examining  Table  II.  The  most  striking  effect  is  on  the  3500°K  point 
with  CO2:  the  revised  data  show  a  substantial  (592‘’k)  Increase  in  adiabatic 
saturation  temperature  over  free  stream  temperature  whereas  the  older  data 
showed  a  significant  (205°)  cooling  effect.  This  is  due  to  a  reduction  in 
the  enthalpy  of  saturated  system  at  a  given  temperature  with  the  revised 
data,  as  shown  by  the  "hsw"  values  given  in  Table  II.  "B"  values  are  not 
greatly  affected,  inasmuch  as  the  revised  data  show  concentrations 

which  are  roughly  one-third  of  the  WO^  concentration  computed  earlier. 

This  is  shown  in  Table  III. 

Comments  on  the  individual  systems  follow. 


BOF-W 

This  system  is  all  condensed  at  1500  and  2000°K  at  1000  psia 
(WOF^/c  and  i^O^/c)  .  disappears  above  2000°K.  The  polymeric 

tungsten  oxide  vapor  species  are  not  particularly  lnq>ortant.  The  reactivity 
Increases  with  temperature. 

CO-W 

At  1500°K  and  1000  psia  this  system  is  all  condensed,  forming 
WC/c  and  A  gas  phase  appears  between  1500  and  2000°K  and  reactivity 

decreases  to  nearly  zero  at  about  4000°K  at  which  point  it  begins  to  in¬ 
crease  due  to  simple  volatilization  of  tungsten.  These  results  are  of 
course  sensitive  to  the  thermodynamic  data  used  for  WC/c,  which  are  not 
well  known.  No  gas  phase  tungsten  oxide  is  significant. 

CO2-W 

This  system  also  forms  WC/c  and  MO^/c  at  1500°K.  WC/c  disappears 

below  2000°K;  MO^Ic  disappears  below  2500°K.  The  polymeric  tungsten  oxide 

species  are  primary  contributors  to  the  ablation.  (WO.(g)  was  indicated 

^  2 

to  be  the  primary  contributor  in  the  last  Quarterly  Report  .) 
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TABLE  II 


COMPARISON  OF  RESULTS  OF  ABLATION  COMPUTATIONS  USING 
"OLD"  AND  REVISED  JANAF  THERMODYNAMIC  DATA  FOR  TUNGSTEN  OXIDES 

1000  psla 


Free  Stream 

Wall  Temperature, 

B  Parameter 

Temp.,  K 

Revised 

"Old" 

Revised 

System  W-BOF 

3000 

2962 

2969 

.0307 

.0292 

3500 

3412 

3443 

.0631 

.0526 

4000 

3901 

3947 

.1275 

.0943 

4500 

4459 

4478 

.2257 

.1563 

5000 

5021 

5004 

.3098 

.2234 

System  W-CO 

2000 

2057 

2057 

.0773 

.0773 

3500 

3501 

3501 

.0009 

.0009 

4000 

4000 

4000 

.00004 

.00004 

4500 

4498 

4498 

.0008 

.0008 

5000 

4981 

4981 

.0065 

.0065 

System  W-CO2 

3000 

2843 

2804 

1.008 

1.148 

3500 

3295 

4092 

1.153 

1.240 

4000 

4611 

4659 

1 .227 

1.141 

4500 

5232 

4984 

1.241 

0.984 

5000 

5624 

5312 

0.911 

0.814 
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TABLE  III 

COMPARISON  OF  EQUILIBRIUM  GAS  MOLE  FRACTIONS 
CO2  plus  W,  1000  psla 

"Old"  and  "Revised"  as  in  Table  II 


Species 

3000° 

Old  Revised 

3500° 

Old  Revised 

4000° 

Old  Revised 

CO 

.7222 

.7781 

.6704 

.8058 

.6947 

.8032 

CM 

0 

0 

.1839 

.1401 

.1051 

.1035 

.0644 

.0829 

0 

.0001 

-- 

.0017 

.0014 

.0099 

.01104 

wo 

— 

-- 

— 

-- 

.0004 

.0004 

.0001 

.0001 

.0009 

.0006 

.0031 

.0038 

WO2 

.0001 

-- 

.0011 

.0008 

.0045 

.0056 

WO3 

.0912 

.0005 

.2203 

.0037 

.2231 

.0145 

W3O9 

.00238 

.0656 

.0004 

.0730 

-- 

.0709 

V12 

-- 

.0155 

— 

.0111 

-- 

.0076 

WO3/C* 

.169 

X** 

1.190 

1.137 

1.334 

1.120 

1.326 

1.143 

hsw*** 

-40.16 

-41.52 

-15.02 

-30.77 

-4.15 

-18.67 

*  Moles/100  grams  of  system  on  an  excess  tungsten-free  basis 

**  Moles  gas  on  same  basis 

***  In  kcal/100  grams  on  same  basis  +  ^298^1-^ 
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2.3  REFRACTORY  OXIDES  FOR  USE  IN  OXYGEN>CONTAINING  ATMOSPHERES 

Equilibrium  calculations  have  been  made  during  the  course  of  this 
work  on  a  variety  of  refractory  materials.  Including  graphite  and  tungsten 
as  well  as  certain  carbides  and  Jjorldes.  All  of  these  materials  have  been 
found  to  be  subject  to  fairly  severe  attack  In  oxygen -containing  atmospheres 
(such  as  CO2  and  H^O) ,  a  fact  which  leads  Immediately  to  consideration  of 
refractory  oxides  as  materials  for  such  applications.  Unfortunately,  the 
oxides  generally  melt  at  much  lower  temperatures  than  do  the  other  materials 
of  Interest  so  that  they  are  not  obviously  attractive  for  high  temperature 
application.  This  Is  shown  by  the  figures  In  the  following  list: 


Compound 

TaC 

HfC 

W 

C 

•  HfN 

HfB2 
Th02 
Hf02 
MgO 
BeO 


Melting  Point,  °K 

4100 

4160 

3650 

^4100  (sublimes  at  1  atm) 
3580 
3520 
3323 
3170 
3073 
2823 


The  highest  melting  oxide  Is  Th02,  which  melts  at  a  temperature 
of  3323°K  or  5521°F,  a  temperature  far  below  that  expected  with  many  ad¬ 
vanced  propellant  systems  of  Interest.  There  did,  however,  appear  to  be 
some  possibility  that  these  materials  would  be  self-cooled  by  their  own 
vaporization  processes,  which  If  appreciable,  could  make  their  usefulness 
somewhat  greater  than  would  be  Immediately  apparent.  A  brief  study  was  thus 
made  of  the  cooling  effect  Involved  with  CO^  as  a  nominal  gas  using  the 
refractory  oxides  MgO,  Hf02  and  Th02.  MgO  was  studied  because  It  Is,  In 
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the  region  of  Interest,  the  most  volatile  of  the  high  melting  oxides; 

HfO^  because  It  Is  the  last  volatile^,  and  of  course  Th02  because  It  Is 
the  highest  melting  oxide. 

The  procedure  used  was  the  same  as  that  employed  In  other  cases 

studied,  except  that  attention  was  devoted  to  the  cooling  effect  during 

adiabatic  saturation  as  a  function  of  pressure.  Thermodynamic  data  for  the 

MgO  system  were  based  on  the  JANAF  tables  to  the  melting  point  (3073°K)^ 

at  which  a  heat  of  fusion  of  18.5  kcal/mole  was  assigned  and  above  which  a 

constant  heat  capacity  of  14.0  cal/mole*°K  was  used.  The  data  for  HfO. 

2  ^ 

and  Its  compounds  was  that  reported  previously  In  this  work  .  These  two 

1  2 

cases  were  computed  using  the  Aeronutronlc  machine  ablation  program  ’  ; 

the  ThO.  case  was  considered  only  by  relatively  crude  hand  computation 

^  8 
methods  Inasmuch  as  the  thermodynamic  data  available  for  Th02  have  not 

been  worked  up  Into  a  form  suitable  for  Inclusion  In  our  machine  program. 

Equilibrium  computation  results  for  MgO  and  Hf02  are  shown  for 

the  three  cases  studied  In  Figure  9  through  13.  Examination  of  the  "A"  and 

"B"  plots  reveals  that  relatively  little  cooling  effect  can  be  anticipated 

In  Che  vicinity  of  the  melting  point,  and  that  where  cooling  Is  substantial, 

saturation  values  are  high.  The  cooling  effect  In  the  CO2  case  Is  made 

more  clearly  evident  In  Figure  13;  a  similar  plot  for  the  Hf02  case  would 

show  substantially  less  deviation  from  Che  45°  line.  (At  3500*^K  and  10  psla, 

the  cooling  effect  with  Hf02  Is  only  5°K.) 

The  self-cooling  effect  due  to  volatilization  of  the  oxides  Is 

summarized  along  with  "ff*  values,  at  the  melting  point,  for  the  three  cases 

studied  In  Table  IV.  The  cooling  effect,  which  Is  shown  In  the  final  column, 

denotes  the  difference  between  the  free  stream  stagnation  temperature  and 

the  wall  temperature  at  the  melting  point  (under,  of  course,  the  assumptions 

used  In  this  work).  Obviously,  the  cooling  effect  Is  greater  at  lower 

pressure  but  In  no  case  Is  It  very  great.  With  MgO,  for  example,  the  free 

stream  stagnation  temperature  can  be  only  approximately  270°  above  the 
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FIGURE  10.  CHEMICAL  CORROSION  OF  HfO^  BY  CO^  AT  THREE  PRESSURES 


FIGURE  11.  SPECIES  PRESENT  UNDER  EQUILIBRIUM  CONDITIONS  AT  THREE 
PRESSURES  FOR  CO  AND  MgO 
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FIGURi  12.  SPECIES  PRESENT  UNDER  EQUILIBRIUM  CONDITIONS  AT  THREE 
PRESSURES  FOR  C0„  AND  HfO 


WALL  TEMPERATURE  (°K) 


AERONUTRONIC  DIVISION 


FIGURE  13.  ADIABATIC  WALL  TEMPERATURE  VS.  FREE  STREAM  TEMPERATURE  FOR 
MgO  IN  CO^  at  three  PRESSURES 
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TABLE  IV 

WALL  COOLING  EFFECTS  WITH  MgO,  Hf02  AND  Th02  AT  THEIR 


MELTING  POINTS 

IN  A  CO2  ATMOSPHERE  AT  VARIOUS 

PRESSURES 

Oxide 

Melting  Point, 

K 

Pressure, 

psla 

••Qii  grams 
ftrams 

Cooling  Effect. 

Or 

MgO 

3073 

1000 

4.1  X  10‘^ 

^3 

100 

6.1  X  10-3 

a^18 

10 

1.4  X  10"^ 

HfO^ 

3170 

1000 

/-O 

<**0 

100 

9.0  X  10'^ 

10 

1.3  X  10"^ 

-/3 

Th02 

3323 

1000 

<^2  X  lO’^ 

^  1 

100 

^3  X  lO'^ 

5 

10 

a»4  X  10‘^ 

'^60 
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3073°K  melting  point  of  MgO  or  3340°K,  even  at  10  psla,  before  melting 
would  be  expected.  With  ThO^  the  maximum  free  stream  temperature  could 
be  only  about  3380°K  even  at  10  psla  before  melting  would  occur;  the 
saturation  parameter  and  thus  the  regression  rate  would  be  smaller  than 
with  MgO.  With  Hf02>  as  expected,  the  cooling  effect  Is  very  slight. 

In  no  case  Is  the  cooling  effect  of  any  real  significance  at  pressures 
of  Interest  In  most  rocket  combustion  systems. 

Based  on  the  foregoing  limited  analysis,  the  oxides  do  not 
appear  attractive  for  high  temperature  oxygen-containing  propellant  systems. 
This  of  course  does  not  preclude  their  consideration  for  all  applications. 

2,4  STABILITY  OF  CARBIDES  IN  NITROGEN 

As  refractory  metal  carbides  are  among  the  highest  melting 
materials  known,  they  are  of  considerable  Interest  as  possible  nozzle 
materials  and  have  been  under  Investigation  In  this  program  In  terms  of 
their  chemical  stability.  These  studies  have  been  limited  In  the  case  of 
the  carbides  reacting  with  nitrogen  In  that  below  certain  temperatures  at 
given  pressures  nitrogen  reacts  quantitatively  with  these  carbides  to  form 
a  condensed  system  consisting  of  metal  nitrides  and  carbon;  the  program 
gives  relatively  little  worthwhile  Information  under  these  conditions.  The 
nitrides  are,  however,  lower  melting  than  the  corresponding  carbides,  as 
shown  In  the  following  tabulation: 

Carbide  Melting  Point, °K  Nitride  Melting  Point, 


HfC 

4160 

HfN 

3580 

TaC 

4100 

TaN 

3363 

ZrC 

3765 

ZrN 

3225 

TIC 

3410 

TIN 

3200 

so  that  It  becomes  Important  to  know  under  what  conditions  carbides  are 

thermodynamically  unstable  In  nitrogen.  (Tungsten  apparently  forms  no 
2 

nitrides  of  Interest  .)  A  brief  study  carried  out  In  this  direction  ensues. 


-27- 


AERONUTRONIC  DIVISION 


The  pressure  of  nitrogen  In  equilibrium  with  the  three  condensed 

1  2 

phases,  C,  nitride,  and  carbide,  was  computed  from  previously  reported  ’ 
and  JANAF  thermodynamic  data,  based  on  the  equation: 

MeC/c  +  ^N2(g)  3^  C/c  +  MeN/c  (1) 

where 

Me  -  Tl,  Zr,  Hf  or  Ta. 

At  the  nitrogen  pressures  and  temperatures  involved  in  the  study, 
the  extent  of  the  nitrogen  dissociation  reaction 

N2(g)  ^  2N(g)  (2) 

and  the  carbide  and  nitride  dissociation  reactions 


MeM/c  ^ 

Me(g)  +  ^2^8) 

(3) 

MeC/c  jdP 

Me(g)  +  C/c 

(4) 

was  found  to  be  negligible  so  that  it  was  necessary  to  consider  only  the 
primary  reaction  (Reaction  1).  Hand  computational  methods  were  thus  used. 

The  pressure  results  are  plotted  against  reciprocal  temperatures 
in  Figure  14.  It  is  evident  from  the  plot  that  at  a  given  temperature  TaC 
requires  the  greatest  nitrogen  pressure  to  convert  it  to  the  carbide  and  HfC 
requires  the  least.  The  relative  stabilities  of  the  carbides  in  nitrogen 
can  thus  be  said  to  decrease  in  the  order  TaC,  TiC,  ZrC  and  HfC.  Even  HfC, 
the  least  stable  of  these,  however,  requires  a  pressure  of  greater  than 
1000  psla  at  2780°K  to  convert  the  carbide  to  the  lower  melting  nitride;  at 
higher  temperatures  (nearer  the  melting  point)  a  pressure  greater  than 
1000  psla  would  be  required.  If  nitride  formation  is  to  occur  in  a  carbide 
nozzle  under  neutral  conditions  in  a  nitrogen  atmosphere,  it  would  form 
thermodynamically  only  during  the  transient  heat  up  period  and  would  not  be 
of  concern  at  the  higher  temperature  conditions  normally  of  Interest.  These 


-28- 


AERONUTRONIC  DIVISION 


(SOOO^K)  (2500“k)  (1000  °K) 

RECIPROCAL  TEMPERATURE,  (°k'^)  ROOlll 

FIGURE  14.  VAPOR  PRESSURES  OF  REFRACTORY  NITRIDES  IN  EQUILIBRIUM  WITH 
BOTH  CARBON  AND  THE  CARBIDE 
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computations  are  of  course  valid  only  for  the  unique  conditions  considered. 
Should  an  additional  reactant  be  present  so  that  solid  carbon  Is  not  present, 
nitrides  may  form  with  specified  nitrogen  pressures  at  considerably  higher 
temperature  conditions. 

For  the  sake  of  completeness,  the  computations  given  In  previous 
reports  for  the  various  carbide -nitrogen  cases  can  be  extended  Into  the 
all -condensed  regions  In  terms  of  the  "A"  saturation  parameter,  using  the 
Information  presented  In  Figure  14  .  These  computations  were  made  for 
1000  psla  nitrogen  pressure  so  that  at  temperatures  below  the  corresponding 
temperature  In  Figure  14  the  carbide  would  be  quantitatively  converted  to 
nitride  and  the  "A"  value  would  be  simply: 

^  _  grams  metal  carbide  _  molecular  weight  of  the  carbide 
gram  nitrogen  14.007 


In  the  case  of  TaC  for  example  at  1000  psla  below  1850  K,  the  A  value 
would  be  simply 

A  .  ^92.96  _ 

^  14.007 


At  temperatures  above  1850°K  at  1000  psla,  the  A  value  would  be  zero  until 
temperatures  at  which  CN  gaseous  species  become  significant.  The  corres¬ 
ponding  temperature  at  1000  psla  for  the  other  systems  would  be  about 
2440°K  for  TIC,  2630°K  for  ZrC,  and  2780°K  for  HfC.  Figure  14  is  useful 
for  other  systems  and  other  pressures  In  a  similar  way.  For  example,  at 
1500*^K  the  equilibrium  nitrogen  pressure  with  TaC  Is  5.5  atmospheres  or 
81  psla.  Below  this  pressure  there  is  no  tendency  for  reaction  to  occur; 
above  this  pressure  the  A  value  would  be  13.8. 

Values  for  the  adiabatic  saturation  parameters  can  not  be  readily 
determined  using  this  simple  approach. 


*  The  equilibrium  composition  plots  In  these  cases  actually  start  at  tempera¬ 
tures  considerably  above  the  all-condensed  region.  Inasmuch  as  finite  quantities 
( ^10~^mole/100  grams)  of  vapor  species  formed  from  the  reactant  solid  must 
normally  be  present  in  the  Aeronutronlc  computer  program  for  computation  to 
proceed.  Thus  where  the  gas  is  1007.  N2  and  the  solid  is  simply  TaC  the  ablation 
program  rejects  the  problem  and  goes  to  a  higher  tenmerature.  In  the  TaC  case  at 
1000  psla,  the  first  machine  results  are  at  3000°K  where  CN(g)  Is  present  at  a 
mole  number  of  5x10"^;  Figure  14  however  shows  that  conversion  to  all  condensed 
occurs  only  below  1850°K. 
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SECTION  3 


PREDICTION  OF  ISOTHERMAL  SATURATION  NUMBERS 
FOR  MIXTURES  FROM  DATA  ON  "PURE"  COMPONENTS , 
WITH  ALLOWANCE  FOR  TOTAL  PRESSURE  EFFECTS 


3 . 1  INTRODUCTION 

In  order  to  extend  the  utility  of  isothermal  and  adiabatic  sat¬ 
uration  numbers  being  determinedunder  this  contract  for  nozzle  materials 
corroded  by  "pure"  combustion  system  components  (CO,  H2O,  BOF,  etc.),  a  ten¬ 
tative  approach  has  been  developed  for  predicting  from  these  numbers  the 
corrosivities  of  complex  mixtures  of  combustion  products  and  at  total 
pressures  other  than  those  reported.  Initially,  the  method  has  been  con¬ 
fined  to  the  estimation  of  Isothermal  saturation  numbers  "A"  for  such  mix* 
tures.  Whereas  the  prediction  of  adiabatic  saturation  numbers  "B"  for 
mixtures  based  on  data  for  the  pure  components  appears  next  to  impossible, 
such  numbers  for  complex  systems  might  be  correlated  empirically,  knowing 
the  general  relationship  between  "A"  and  "B"  for  the  individual  components. 
Accordingly,  corrosion  rates  might  be  estimated  for  actual  propulsion  ap¬ 
plications,  assuming  the  adiabatic  saturation  model. 

The  method  thus  developed  is  not  intended  to  replace  the  rigorous 
calculation  of  saturation  numbers,  using  the  ablation  equilibrium  computer 
program.  Rather,  it  is  designed  to  permit  rapid  estimation  of  "A"  values 
for  cases  where  a  computer  or  suitable  program  is  not  available,  or  in 
which  rigorous  computation  would  require  excessive  use  of  machine  time. 
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The  basis  of  the  approach  first  will  be  outlined,  Including  the 
preliminary  "additivity"  assun^tlon,  and  the  more  elaborate  model  subse¬ 
quently  developed.  Mathematical  treatments  then  will  be  given  for  four 
typical  cases,  and  possible  extensions  to  other  situations  will  be  In¬ 
dicated.  Finally,  two  examples  will  be  presented,  showing  actual  com¬ 
parisons  between  predicted  and  rigorously  computed  Isothermal  saturation 
numbers . 

3.2  THEORETICAL  BASIS 

A  simple  assumption  which  might  be  used  to  estimate  Isothermal 
saturation  numbers  "A"  for  complex  systems  Is  the  rule  of  "additivity".  In 
which  the  overall  "A"  would  be  taken  as  a  linear,  weighted  average  of  the 
"A's"  of  each  of  the  Individual  combustion  species.  In  applying  this  ap¬ 
proach,  the  "A"  for  each  component  would  be  taken  at  the  total  system 
pressure,  rather  than  at  Its  own  partial  pressure,  since  otherwise,  "A" 
values  of  Infinity  might  be  determined  for  some  trace  species.  Even  with 
this  technique,  an  "A"  value  for  a  certain  species  may  be  undefined.  If  a 
question  exists  as  to  whether  the  species  yields  a  gaseous  phase  In  equi¬ 
librium  with  the  ablator  at  the  specified  temperature  and  pressure. 

As  will  be  shown,  the  "additivity"  model  by  Itself  yields  correct 
predictions  of  Isothermal  saturation  numbers  only  for  a  very  limited  number 
of  systems.  Two  obvious  factors  tending  to  Invalidate  this  assumption  are: 
(1)  partial  pressure  effects  on  equilibria  Involving  a  change  In  gaseous 
moles,  and  (2)  chemical  Interaction  between  two  corrodant  species.  The 
method  herein  developed  Is  concerned  primarily  with  applying  "corrections" 
based  upon  factor  (1)  for  the  case  of  dilution  of  the  attacking  species  by 
an  essentially  Inert  carrier  gas.  By  some  modification,  the  approach  can 
be  expanded  to  Include  dilution  by  an  Independently  corrosive  agent,  or, 
under  factor  (2), even  dilution  by  a  second  species  entering  Into  the  same 
equilibrium  reaction.  The  only  situation  not  amenable  to  direct  analysis 
Is  that  In  which  a  new  ablation  species  appears  upon  going  to  the  complex 
system. 
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3.3  A  CATALOG  OF  ABLATION  REACTIONS 

In  order  to  characterize  ablation  reactions  for  possible 
pressure  effects,  four  main  groupings  were  defined,  as  shown  in  Table  V, 
according  to  whether  the  net  mole  change  was  zero,  positive,  negative, 
or  independent  of  the  gaseous  species.  Various  ablation  reactions,  in- 
voling  corrosion  of  C/c,  TiC/c,  ZrC/c,  ZrB^/c  and  W/c  nozzle  materials 
have  been  classified  according  to  this  scheme;  the  results  are  presented 
in  Table  VI.  As  an  approximation,  any  reaction  involving  a  mole  change 
of  less  than  50  percent  (2-^3,  or  3^2),  has  been  considered  as  having 
zero  change  in  gaseous  moles. 


TABLE  V 

CLASSIFICATION  OF  ABLATION  REACTIONS 


Gas -Phase 

Case 

Description 

Type 

Mole  Reaction 

I 

No  change  in  gaseous  moles 

I 

n-i^n 

II 

Increase  in  gaseous  moles 

11 

II' 

II" 

II"  ' 

Il(4) 

1- >2 

2- »-3 
1-^3 

1 — ^4 

and  higher 

0-+n 

III 

Decrease  in  gaseous  moles 

III 

2->l 

III' 

3— ^2 

III' ' 

n-e*0 

IV 

Inert  gas  feed;  simple  de- 

IV 

0->n 

composition  or  vaporization 
of  ablator 
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TABLE  VI:  A  CATALOG  OF  ABLATION  REACTIONS  (Pressure  •  1000  psla)  (page  2) 
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3.4  ANALYTICAL  SOLUTIONS  FOR  TYPICAL  CASES  AND  COMMENTS 

3.4.1  MATHEMATICAL  DEVELOPbCNTS 

Mathematical  treatments  are  presented  In  the  following  (Case  I, 
Case  II,  Case  III  and  Case  IV)  to  show  how  Isothermal  saturation  numbers 
for  typical  ablation  reactions  with  Inert  diluents  are  determined  for  each 
of  the  four  main  classes. 


Case  I;  No  change  In  gaseous  moles 

Reaction:  X/c  +  Y  XY 

Assume:  X/c  Is  essentially  non-volatlle. 

Gaseous  component  Z  Is  Inert. 

a.  For  pure  Y  In  feed: 

■f-'f  "  ” 


b.  For  mixture  of  Y  and  Z  In  feed: 


(1) 


(2) 

(3) 

(4) 


molt  »  ■f  t^elts  Y  *  molts  Z.'j  ^5) 

•  0-  :f)  7T  <8> 
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( t  )  f  Mx 

fMy 


or,  In  other  words 


fMy 


-). y 

'  V  •t'oiat 


ih  -fted  ) 
fted  f 


showing  the  linear  relation  for  additivity  when  no  mole  change  Is  Involved. 


c.  For  mixture  of  Y  and  XY  In  feed: 


Assume  that 


\j  nnole  fe^d  ~  ^  ol^s-  Y  *•  (l-  S)  Xy] 

-  (T^)rr 

^XY  ‘  (TTT-r)^ 

(~frT~r^  >  O-f-),  A  tfv  ^  y  ^;ry  ' 

Ay  *  ^  ~  (15) 

My  -f-  (l-J-)Mx 


A  a  I  )  ~  (f~f)]f^X 

i^My  -h  0-;i)MxY 


or  In  other  words, 


'Y*  XY 


!  Ja,.  -  \'j/Y,,fthf  Y  ill  \ 

/  A  y  Af^y/J\  Vo-Ya!  f^td)  ^ 

ir  [Ay  ~  ^Xmy)]  f'*'  oi/\tfayr  /etr  tMmn  Ay 

ase  although  there  Is^no  change  In  moles  on  reaction,  the  s 


w^*rl•  ~yr  IAy  ~  ' 

In  this  case  although  there 
additivity  rule  would  fall. 


simple 
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Case  II:  Increase  In  gaseous  moles 

Reaction;  X/c  +  WY^^fcW  +  XY  (1) 

Assume:  X/c  Is  essentially  non*volatlle. 

Gaseous  component  Z  Is  Inert, 
a.  For  pure  WY  In  feed: 


^wy  “ 

n-  Z  k[-  !  +-Jl  +  ^  ) 

(2) 

^XY  • 

(3) 

II 

X 

> 

/M)(  \ 

(4) 

L  n  -  /<(-  /  +y/ 

V  My  ' 

For  mixture 

of  WY  +  Z  In  feed; 

[l  mol*  f**d  «  /  molts  WY  *  (l""  molts  Zj 

0  -*-^0  ^^(^)]} 


. 
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or 


where 

Suf  A^y  ^ 

Given; 

A^.y.  f-/  *l//  ] _ <W«  I 

<(f^-[-i*{T^^y^-y' 

and 

4'  .  <2'^f~ (>  *y C/*  j  J _  !  My  \ 

4^)-(i)[-(l*^)*Tf(l*f)‘*-tf(^)]}  'Mw/ 

Then;  _ _ 

^WY  f-  /  -^y /  v-(-g-;  7 

X  jj 

Given  >  one  can  plot 

'n  ' 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


vs .  y 
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M  I  ,  ^ 
A 


WY 


*/  (JdjL.  ) 
At  f  ^  o .  0^  ar 

'  A^^Y  A 


WY 


[ 


Ot 


A'  «  7 


tdhge  /  *  ^  /'o  /  ■  /^  >  / 


'  wy 


c.  For  mixtures  of  WY  +  W,  or  WY  +  XY  in  feed: 

(This  case  not  derived.) 


(17) 


(18) 


(19) 


(20) 
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Case  III!  Decrease  In  gaseous  moles. 

Reaction:  X/c  +  2Y  XY^ 

Assume:  X/c  Is  essentially  non-volatile. 

Gaseous  component  Z  is  inert. 


a.  For  pure  Y  feed: 


(1) 


(2) 


TV  - 


(3) 


and 


(4) 


b.  For  mixture  of  Y  and  Z  In  feed: 

/"/  me/«  “  /  mol^s  Y  +  0“^)  rr>ol(fS  7L  J  (5) 

/■y  = 

Ay, '  ^  m 

■Pz  «  (/ -/J p  "  <*> 


and 


'y-fZ 


^/[l*£f&-j<)Kn]-^/l*4A£-/)Kn  If  Mx  I 

[[I  *4(a-f)KT\]-^ll  *4H^-f)Kn  *  0-f)M2,J 


(9) 
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where 


Given  Kt,  one  con  plot  ^X-  vs. 

Ay 


Ai-  /■*  /^  ^  a  / 


(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
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Case  IV:  Inert  gas  feed;  simple  decompoflltton  of  ablator. 

Reaction:  X/c  X  (1) 

Assume  and  Z2  are  Inert. 

a.  Feed  composition:  1  mole  feed  >  f  moles  and  (1-f)  moles  Z^.  (2) 

•>  H  (3) 

.  f(n-K)  (4) 


But 


^2,  -  ' 
A  s  /  ^tL. .  \  M 


Thus 


•Z, 

MiL 

AZjAz^ 


Of- 


(6) 

(7) 

(8) 

(9) 


(10) 


The  additivity  rule  applies  In  this  case,  even  though  there  Is  a  change  In 
gaseous  moles. 


b.  Feed  composition:  X  +  Inert. 

(This  case  not  derived.) 
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3.4.2  COMMENTS 

It  is  seen  that  the  simple  additivity  rule  applies  only  for  re¬ 
actions  Involving  zero  mole  change,  or  for  ablator  decomposition  reactions 
with  Inert  diluents.  For  reactions  Involving  an  Increase  In  gaseous  moles, 
equilibrium  ablation,  as  would  be  expected.  Is  enhanced  by  either  a  drop 
In  total  pressure  or  dilution  with  Inerts;  exactly  the  reverse  Is  true  for 
reactions  In  which  net  gaseous  moles  disappear.  As  part  of  Case  I,  It  Is 
also  shown  that  dilution  with  an  ablation  product  species  reverses  the 
corrosion  equilibria,  yielding  lower  Isothermal  saturation  numbers  than 
dilution  with  a  comparable  amount  of  Inerts.  The  same  should  hold  true 
for  the  other  three  cases  as  well. 

Within  certain  limitations,  the  method  may  be  extended  to  cover 
other  situations  as  well.  For  reactions  having  mole  changes  such  as  2-*  3, 
1—^3,  or  1-^4,  rigorous  mathematical  treatment  becomes  unwieldy,  since  a 
third -or -higher  order  algebraic  equation  must  be  solved.  It  Is  obvious, 
from  extension  of  Case  I,  that  the  "addltlvlty"rule  still  applies  to  two  or 
more  Independent  ablation  reactions.  If  the  mole  change  Is  zero  for  each. 

If  the  mole  changes  are  not  zero,  the  problem  of  two  Independent  reactions 
Is  much  more  complex.  Fortunately,  one  ablation  equilibrium  usually  pre¬ 
dominates  over  the  others  In  mixtures  of  practical  Interest. 

3.5  EXAMPLES 

Two  examples  are  provided  to  Illustrate  the  application  of  the 
method,  and  to  point  out  some  calculations!  details  not  apparent  from  the 
preceding  discussion.  The  first  system  involves  BOF  as  the  major  corrodant, 
while  In  the  second  case,  H2  la  the  principal  reacting  species.  In  both 
examples  the  ablator  Is  graphite  (C/c). 
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Example  1:  ■  3865°K,  ■  500  psla,  C-H-N-O-F-B  propellant. 

The  computer-derived  equilibrium  composition  of  the  free  stream 
at  the  given  temperature  and  pressure  was  recalculated  In  terms  of  the 
minimum  number  of  stable  species  required  to  make  up  the  system  at  these 
conditions.  Since  the  Isothermal  saturation  number  "A"  Is  fixed  by  overall 
equilibrium  rather  than  by  reaction  rates,  this  simplification  appears  a 
valid  means  of  eliminating  minor  or  uncertain  species  from  consideration. 
For  the  system  under  consideration,  the  composition  was  most  conveniently 
expressed  in  terms  of  the  following:  BF,  BOF,  CO,  HF,  N^. 


Composition _  ■^1000 » 


Gm-moles 

100  sms 

Gms 

100  sms 

Isbthermal  Saturation 
at  1000  psla 

Niffliber 

BF 

0.050 

1.49 

•k 

0 

BOF 

0.872 

39.93 

0.27 

CO 

1.286 

36.02 

0.0035 

HF 

0.021 

0.42 

0.17 

«2 

0.088 

0.18 

5.0 

0.784 

21.96 

0.072 

* 

Estimated  from 

"A's"  of  BFj  and 

BF2. 

On  preliminary  examination,  the 

major 

corrodant  species 

is  seen 

!  BOF 

,  based  on  both 

its  "A"  value  at 

1000  psla,  and  upon  its 

large 

weight  fraction.  (H2  has  a  much  higher  value  of  "A",  but  based  upon  its 
very  small  weight  fraction,  should  not  be  a  significant  reactant.)  The 
ablation  reaction  of  BOF  with  graphite  is: 


C/c  +  BOF 


CO  +  BF 


-51- 


ACRONUTRONIC  DIVISION 


This  equation  fits  the  form  of  Case  II: 

K 

X/c  +  WY  Sit  W  +  XY  (Eq.  II-l) 


In  which  X  Is  C,  W  Is  BF,  and  Y  Is  0.  The  equilibrium  constant,  may 

be  obtained  from  the  JANAF  thermochemical  data  tables,  or,  as  an  alternative, 
evaluated  from  the  computed  "A"  value  for  BOF  at  3865°K.  '  In  order  to  allow 
for  the  effects  of  minor  side  reactions,  the  latter  Is  preferable: 


0.27  = 


-1  + 


(iLOi) 

'■45.82'' 


(Eq.  II-4) 


from  which  =  oo  ;  l.e.,  the  ablation  reaction  goes  essentially  to 
completion. 

Knowing  that  the  "A"  value  for  BOF  Is  determined  by  stoichiometry 
rather  than  by  equilibrium,  the  approach  must  be  modified  as  follows: 

(1)  The  value  A  ■  0.27  for  BOF  at  SEbS^K  can  be  assumed 
constant.  Irrespective  of  the  total  pressure  or  the 
dllut Ion. 

(2)  The  Influence  of  the  large  CO  concentration  on  the 
BOF  ablation  reaction  can  be  neglected. 

(3)  Since  the  BOF  stoichiometry  Is  fixed,  some  other 
species  can  be  considered  as  a  major  corrodant. 

The  next  most  corrosive  species  Is  seen  to  be  N^,  by  virtue  of  the 


reaction: 


2C/c  +  N2 


2CN 


The  equation  also  fits  the  form  of  Case  II: 

K 

X/c+WY^  W  +  XY  (Eq.  II-l) 


In  which  X  Is  2C,  W  Is  CN,  and  Y  Is  C’l-N. 
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As  before,  the  equilibrium  constant,  "K2",  may  be  evaluated  from 
the  computed  value  of  "A"  at  3865°K  and  1000  psla: 


0.072  =: 


-1  +  y  1  +  (^) 

K2 


.24.02. 

^28.02^ 


(Eq.  II-4) 


from  which  K2  "r  0.49.  At  500  psla  pressure  and  3865°K,  the  "A”  value  for 
pure  N2  attacking  graphite  would  be: 


A  — 
*500  - 


-1^ 


(*Thls  computation  shows  how  the  method  can  be  used  to 
compute  total  pressure  effects  but  Is  not  otherwise 
used  In  the  material  which  follows.) 


In  the  given  mixture,  the  mole  fraction  of  is  0.252.  Accordingly,  the 
"effective"  Isothermal  saturation  number  of  In  the  mixture  may  be  derived: 

\  f  -1.252  +  /(1. 252)^  (4)  (0.252) 

,  w  _ L _ * 

0.252  |(5^)  -  1-1^-1.252  +  /(1. 252)^  +  0 

(Eq.  11-15) 

Neglecting  effects  of  pressure  and  dilution  upon  the  ablation  equi¬ 
libria  for  the  minor  corrodent  species,  the  overall  "A"  value  for  the  com¬ 
bustion  mixture  at  3865*^K  and  500  psla  can  be  estimated: 


0(0.252)  ( 


oArJ 


(24^) 

'•28.02'' 


0.217 


A  =  [_(1.49)(0)  +  (39.93)(0.27)  +  (36.02) (0.0035)  +  (0.42)(0.17) 

+  (0.18)(5.0)  +  (21.96)(0.217)J  -r  100  -  0.166 

This  compares  to  an  "A"  value  of  0.168  determined  rigorously,  using  an 
IBM  7090  computer. 
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Example  2;  -  2793°K,  ■  1000  psla,  B-N-H  propellant 

As  in  Example  1,  the  computed  equilibrium  composition  was  re¬ 
calculated  in  terms  of  the  minimum  number  of  stable  species.  For  the 
system  under  consideration,  the  conq>ositlon  was  expressed  in  terms  of 
N2,  and  BN/c: 


_ Composition _ 

Gm-moles  Gms 

100  gms  100  gns 


Aiooo* 

Isothermal  Saturation  Number 
_ at  1000  psia _ 


H2  6,609 
N2  0.165 
BN/c  3.305 


13.32 

4.62 

82.06 


0.44 

0.0009 


(^Assumed  non-reactive.) 

Since  BN/c  was  considered  inert,  calculations  were  based  only  on  the  gas  phase. 

The  major  corrodant  species  is  seen  to  be  H2.  At  this  temperature 
and  pressure,  the  primary  ablation  reaction  of  H2  with  graphite  is: 

K 

2C/c  +  H2  C2H2 


This  equation  fits  the  form  of  Case  1: 

K 

X/c  +  XY  (Eq,  I-l) 

in  which  X  is  2C  and  Y  is  H2. 

Because  the  reaction  involves  no  change- in  gaseous  moles,  dilution 
effects  can  be  neglected. 

The  isothermal  saturation  number  of  only  the  gas  phase  may  be 

estimated: 


.  -  (13.32)(0.44)  +  (4.62) (0.0009) 

gas  “  (13.32)  +  (4.62) 


0.327 
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Based  upon  the  total  combustion  mixture,  the  Isothermal  saturation  number 
would  be  predicted  as: 


(13.321 


(0.0009) 


f82.06) 


0.059 


This  compares  to  an  "A"  value  of  0.079  determined  rigorously,  using  an 
IBM  7090  computer.  In  this  Instance,  the  method  does  not  give  as  good  an 
estimate  of  "A",  as  desired. 

By  Inspection  of  the  computed  equlllbrltim  composition  of  the  com¬ 
bustion  mixture  after  saturation  with  graphite,  the  reason  for  the  discrepancy 
between  predicted  and  computed  "A"  values  was  uncovered.  A  major  product  of 
the  ablation  process  was  found  to  be  gaseous  HCN,  which  obviously  was  not 
taken  Into  account  In  the  Individual  "A"  values  of  and  with  C/c  con¬ 
sidered  separately.  Owing  to  a  three-way  chemical  Interaction  between  two 
combustion  species  and  Che  ablator,  the  prediction  approach  breaks  down 
In  this  case. 
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SECTION  4 

PLANS  FOR  FUTURE  WORK 


It  Is  contemplated  that  future  effort  under  this  contract  will 
be  devoted  to  the  following: 

(1)  Detailed  study  of  HfN,  the  highest  melting  nitride, 

In  equilibrium  with  various  "pure"  gases. 

(2)  Interpretation  of  computer  results  to  compare  suita¬ 
bility  of  various  materials  In  various  atmospheres, 
considering  melting  points  and  density  effects.  The 

degree  of  cooling  [  8tream-^wall>/® J 

unit  of  wall  material  evaporated  or  In  other  fashion 

will  be  considered. 

(3)  Comparison  of  various  materials  for  real  propellants. 

(4)  Comparison  of  theoretical  corrosion  rate  predictions 

against  measured  values  for  real  propellants  or  pure 
gases.  * 

Certain  other  areas  such  as  solid  (or  liquid) -solid  reactions, 
effects  of  uncertainties  In  thermodynamic  data,  further  mixture  and  pressure 
effects,  and  a  general  study  of  other  types  of  materials  will  be  carried  out 
as  time  permits. 
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APPENDIX 

ERRATA  IN  PREVIOUS  REPORTS 


A  number  of  significant  errors  or  omissions  In  presentation  have 
been  noted  In  the  three  previous  quarterlies  Issued  under  this  contract, 
as  follows : 

First  Quarterly  (15  September.  1962) 

p,  16  --  should  be  defined  as  moles  of  nozzle  material  as  condensed  phase 
per  100  grams  of  equilibrium  mixture,  Including  condensed  phase, 
p.  17  and  Figures  1-8  --  the  pressure  used  In  these  calculations  was  1000  psla. 

Second  Quarterly  (15  December.  1962)  -  Volume  I 

p.  66  — ^*^£298  ^^5(8)  should  be  listed  as  -384  kcal/mole. 

Third  Quarterly  (15  March.  1963) 

p.  4  --  Brackets  should  surround  the  terms  following  T,  as 
[Ht  -  Hj,,  -  <T-298)J 

p.  6  —  The  symbol  Identification  code  should  read 

o  Southern  Research 
a  Mezakl 
X  K.  K.  Kelley 

p.  33  --  The  title  for  Figure  9  should  have  the  words  "at  1000  psla"  added, 
p.  36  --  The  title  of  the  referenced  publication  should  read 

"Thermodynamic  Properties  of  Inorganic  Substances,  V.  High 
Temperature  Heat  Contents  of  Fifteen  Refractory  Borides." 
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